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Abstract 

We estimate the long distance (LD) contribution to the magnetic part 
of the s — > d'y transition using the Vector Meson Dominance approximation 
(V = p,cu,ipi). We find that this contribution may be significantly larger 
than the short distance (SD) contribution to s —> dy and could possibly 
saturate the present experimental upper bound on the O - —> decay 

rate, ~ 3.7 x 10” 9 eV. For the decay B s —*• B^, which is driven by 

s —» dy as well, we obtain an upper bound on the branching ratio BR(B S —> 
B* d 7) < 3 x 10" 8 from r^ s - r 

Barring the possibility that the Quantum Chromodynamics coefficient 
a 2 (m s ) be much smaller than 1, also implies the approximate re¬ 
lation | W + \^r-- This relation agrees quantitatively with 
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a recent independent estimate of the l.li.s. by Deshpande et ah, confirming 
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that the LD contributions to b —> sj are small. We find that these amount to 
an increase of (4±2)% in the magnitude of the b —> sy transition amplitude, 
relative to the SD contribution alone. 

1. Introduction and Overview 

The investigation of the quark radiative transition b —> sy has been an 
important focus of attention in recent years [ 1 ] both because of experimental 
measurements [2] and because long distance (LD) corrections to the Stan¬ 
dard Model (SM) predictions for the short-distance (SD) contributions are 
estimated to be small [3]. (For exclusive B —> K *7 decays see Ref. [4]). Thus, 
this transition constitutes an excellent laboratory to test the SM or possible 
high energy deviations thereof [5]. It has been pointed out recently [ 6 ] that 
for the c —» wy transition the situation is reversed, with the LD contributions 
dominating over the SD ones by many orders of magnitude. 

In this paper we investigate the analogous quark transition s —> dy and 
two exclusive hadronic processes, —> S _ y and B s —> R*jy, where it plays 

an important role. Throughout this paper we are concerned with the mag¬ 
netic transition only, since the charge-radius one vanishes for real photons. 

The SD contribution to s —> dy has been investigated before [see e.g.: 
7,8,9] and we simply repeat the calculations, using updated values for the 
relevant QCD coefficients. Applying the quark model formalism of Ref. [10] 
we find that the SD s —> dy contribution (by itself) to the decay 
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rate is far below (by a factor of order 600) the present experimental upper 
limit [ 11 ]: 

T(Q~ -> S- 7 ) < 3.7 x 10 ~ 9 eV (90%CL). (1) 


Hadronic LD effects that involve light mesons in loops are estimated to 
be small[ 8 ,12], comparable to the SD contributions. 

On the other hand, by using a Vector Meson Dominance (VMD) approx¬ 
imation for the LD contribution to the s —> dy transition (along the lines 
discussed by Deshpande et al. [3] for b —> sy), we find LD contribution that 
are likely to be significantly larger than the SD ones. In fact, the rate for 
O - —> may not be far from the experimental bound (1), due to this VMD 
contribution. The resulting VMD amplitude is approximately proportional 
to 


a 2 {m s 


2 y 4 (°) 1 3p(°) lfl£(0) 

3 V 


mz 


2 m 2 
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where a, 2 (m s ) is a Quantum Chromodynamics (QCD) coefficient [13] and the 
< 7 y( 0 )’s are the usual vector meson-photon couplings, evaluated at q 2 = 0 . 
Although a direct estimate of a, 2 (m s ) is not reliable because we are well 
into the low energy region where perturbation theory cannot be trusted, we 
can use ‘“smoothness” arguments to extrapolate from the phenomenologi¬ 
cally determined values a 2 (ml) = 0.24 ± 0.04 [3] and a 2 (m^) = 0.55 ± 0.1 
[13] to obtain a, 2 (m s ) > 0(0.5). We then apply the formalism of Ref. [10] 
to get an expression for the fV —> decay rate from our SD+VMD 
s —> d'y amplitude. (Notice that there are no pole contributions to this de¬ 
cay). It turns out that if the above rough estimate a 2 (m s ) > 0(0.5) is correct 
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then the experimental limit ( 1 ) can be satished only if the contribution of 
the V-t; resonances in the parenthesis of eq. (2) cancels, at a level of 30% or 
better accuracy, the p and u meson contributions, which can be reliably ob¬ 
tained from the p and c o leptonic widths [14], The limit (1) then forces the 
approximate relation at q 2 = 0 : 

t3l 

3 V rn% 2 w? p 6 ml 1 J 

which is highly nontrivial, and may be interpreted as a remnant of the badly 
broken SU(4)p symmetry. 

The relation (3) turns out to be very useful for the b —> sy decays. As 
noted in Ref. [3], in the VMD approximation the main LD contributions to 
this decay can be expressed in terms of the l.h.s. of eq. 3. In Ref. [3], the 
sum in the l.h.s. of eq. (3) is estimated by using measured leptonic widths 
of the states and ^ photoproduction data as well as an assumption about 
the higher ^ excitations. We estimate this sum with better accuracy by 
replacing experimental values for g p (0) and g u ( 0) in the r.h.s. of eq. (3) and 
find very good quantitative agreement with the central value obtained in 
Ref. [3]. We thus confirm the main result of Ref. [3] that LD contributions 
to b —> s'j are of order of a few percent. According to our explicit estimate, 
these corrections amount to an increase of (4 ± 2%) in the magnitude of the 
b —> sy transition amplitude, relative to the SD contribution alone. 

Finally, we also apply the SD + VMD approximation for s —> d'y to the 
unusual decay mode B s —> B ^ 7 , where the b quark plays the “spectator” role. 
We point out that this decay (followed by B* d —> B ^ 7 ) has a clear experimen¬ 
tal signature of 2 monochromatic photons of energies ~ 50MeV each. We 
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find, using the limit of eq. ( 1 ), a small but hopefully measurable branching 
ratio BR(B S —> B*tf) < 3 x 1CT 8 . 

2. SD Contribution to the s — > dy Amplitude 

The SD amplitude relevant to the s —> dy transition can be expressed as 

c G 

A S d = - g ^-2 ^ F 2 (/^ 2 )da^[m s -R + md-Ljs-F^ , (4) 

where m a , rn c i are current quark masses and F 2 (/j' 2 ) is a form factor evaluated 
at a low scale // > 0(m s ) which includes (dominant) QCD corrections. Early 
estimates of F 2 (yU 2 ) [7,8,12] were in the approximate range 0.15 — 0.36 [15] 
while we obtain by explicit calculation, using ay(m 2 ) ~ 0.3, a s (fi 2 ) = 0.9 
in the formulas given in Ref. [16], a somewhat smaller value F 2 (yU 2 ) ~ 0.1, 
which will be used below (see also Ref. [9]). 


3. LD Contribution to s —> dy 

To estimate the LD contribution to s —> dy we use the VMD approxima¬ 
tion in analogy to the formalism used in Ref. [3] for b —> sy. As an intermedi¬ 
ate step one defines a transverse amplitude A(s —> dV(q))x (V = ^,p, u in 
this case) and then introduces the V to 7 conversion vertices, setting q 2 = 0. 
Using Gordon decomposition we find that the LD amplitude for the s —> dy 
transition is 


At, n = 




'2 gim 

, 3 , m 2 


1.9p(0) 1^(0)' 
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where we have used V cs V* d ~ —Vu S V* d , a 2 (/i 2 ) is a QCD coefficient the value 
of which is taken from phenomenology in the context of the factorization 
approximation [13], and the gv(q 2 ) factors are defined in the usual way, e.g. 
(^(g)|c7 M c|0) = igp(q 2 )e+(q). We have not included possible contributions 
from the p and u radial excitations ( p',p",... ...) because we think 

that their contribution is much smaller and is already taken into account 
to a significant degree in the SD amplitude (4). The i/j excitations should 
be included however, because they are narrow resonances that are clearly 
distinguished from the cc continuum. Note that due to the hadronic nature 
of the VMD approximation, M s and M d should correspond to “constituent” 
mass parameters. (The use of “constituent quark” spinors in deriving (5) 
should take into account to some extent non-perturbative effects such as chi¬ 


ral symmetry breaking and confinement). In any case, it turns out that only 
the combination which has a similar magnitude for “constituent” or 

“current” s, d quark masses, appears in our applications (see Sects. 4 and 5) 


when the interference between the (presumably) dominant LD contribution 


and the SD contribution is neglected. 


It is difficult to estimate the coefficient a 2 (p 2 ) for p > 0{m s ) appearing 


in eq. (5). However, a smooth extrapolation from the phenomenologically 
obtained values a 2 (m 2 ) ~ 0.24 ±0.04 and a 2 (m 2 ) = 0.55 ±0.1 [3,13] leads to 
a 2 (m 2 ) > 0.5. 

The couplings gp z (m^), g p (m 2 p ), g^m 2 ) are readily obtained from lep- 
tonic decays of these mesons, but their extrapolated values at q 2 = 0 are 
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less trivial, especially for the 7/y states. Photoproduction data seems to in¬ 
dicate that g 2 (0) ~ g 2 (m 2 p ), g 2 ( 0) ~ g 2 (ml,) [17,18]. On the other hand, 
estimates in Ref. [3] using t/> photoproduction data [18-20] give g 2 ( 0) = 
(0.12 ± 0.04)^(777.^,). In Ref. [3] it is also assumed that the same ratio holds 
for the excitations i/j 1 , t//', etc. 

Making use of the above estimates as well as of the lcptonic widths of 
the relevant vector mesons [14] we obtain the numerical values g 2 (0)/rn 2 ~ 
0.047GeV 2 , (? 2 (0)/m 2 ~ 0.038GeV 2 and Ei - — 0.041GeV 2 . The hrst 
two estimates should be accurate to about 10% while the latter must be 
considered only as a rough estimate, with an uncertainty of at least 40%. 
Once we derive the approximate relation (3) we will be able to give a far 
more reliable estimate of V, - . which is consistent with the above central 
value. 

4. Application to the Decay and Consequences 

We use the quark model of Ref. [10] to estimate the rate for the decay 
0“ —> S - 7 , from the SD and LD contributions to the s —> dy quark decay 
amplitude obtained in previous sections. 

For notational convenience, we define the constants v = [RcsKdl — 0-22 
and C V md = (f E; • The relative sign of the SD and 

LD contributions is determined by the theory [3] so that the full amplitude 
for the s —> dy transition can be written as 
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Atot(-s 


d'y) — A sd + A ld 
tGp -dar 


' m s F 2 va 2 C\r M i)M s 

. 8tt 2 + 


Ml - Ml 


V2 

fm<iF 2 va 2 Cy\,mM d \ \ 

Ml - Mj ) ) S ^ ' 
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Following Ref. [10] we then obtain 


r(ir 


- 7 ) = 


aG 2 F fm = 



I# 

87r 2 na2CvMD^s \ 

M? ~ Ml J 


( 6 ) 


+ rn d F 2 - 


87T 2 Va 2 C'vMD M d 

Ml ~ Ml 


where <f is the photon momentum in the fl _ rest frame and the separate SD 
and LD contributions are exhibited explicitly. 

In the absence of LD (VMD) contributions, we would obtain (for m s ~ 
175MeV, rn d ~ lOMeV, F 2 ~ 0.1, see Sect. 2) 


Tsd(o 


"7) ~ 6.4 x 10 2 eV 


( 8 ) 


which is far below the present experimental bound of r exp (0 —> H 7) < 
3.7 x 10~ 9 eV. On the other hand, the large theoretical uncertainty of over 


40% in the value of the sum 


alX 0) 




(see Sect. 3) which appears in Ovmd, 


would allow the LD contribution to saturate this experimental bound. In fact, 


the experimental limit can be used to constrain Ovmd and hence 
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Using typical values M s ~ 0.5GeV, M ( [ ~ 0.35GeV for the constituent quark 
masses and 02 > 0.5 (see Sect. 3), we find 


\C- 


VMD 


|2 y. 4i( Q ) _ 1 ^( Q ) _ 1 ^(°) 
, m3,. 2 m3 6 m 2 , 


< O.OlGeV 2 


(9) 


V * “ " “p 

This constraint would be only slightly different, had we used current quark 
mass parameters instead of M s and M ( j. 

The bound in eq. (9) represents a remarkable cancellation at the 30% level, 
considering that + \ ^$r- — 0.030GeV 2 (see Sect. 3). We presume that 
this effect may stem from the combination of the GIM [21] mechanism and the 
underlying SU(4)^ symmetry, which if exact would give a full cancellation 
(after inclusion of p',pf' a/, u" ... states). The SU(4)^ symmetry is 
known to be badly broken by the large mass of the c quark. However, here 
we are comparing the form factors g^.(q 2 ), g 2 (q 2 ), g 2 (q 2 ) at a common scale 
q 2 = 0, which seems to “restore” this symmetry to some extent. We have 
noticed that if |^(0)| ~ |^(m|)| [17,18], leading through 0 leptonic width 
data [14] to |<?^(0)| ~ 0.24GeV 2 , a completely analogous near cancellation 
occurs for the quantity C0 MD = — \^r + — which is relevant to 

<f> P & 

LD effects in c —> u'y decay [6]. We obtain C0 MD ~ —1.8 x 10” 3 GeV 2 , which 
represents a cancellation at a level better than 10% for which presumably 
the SU(3)i? symmetry is responsible. 

We note that the upper bound (9) on |CVmd| tells us that although the 
LD effects are likely to dominate s —> d'y, they can be at most a factor of 
about 25 larger than the SD contribution in the amplitude. This represents 
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an intermediate situation between the b —> sj decays where the SD contribu¬ 
tion clearly dominates [3,4,22] and the c —» u'f decays where the SD effects 
are completely negligible relative to the LD ones. 

5. Implications for the LD Contribution to b — > sy 

Because — 0.030GeV 2 , eq. (9) implies that the approx¬ 

imate relation given in eq. (3) must hold to an accuracy of order 30%. This 
then independently determines 

x;^ = 0 . 045 ± 0 . 016GeV 2, (io) 

i 

where our uncertainty in the values of g p (0) and < 7 W ( 0 ) has been folded in. 
Notice that this result is in very good agreement with the central value 
(~ 0.041) estimated from t/j photoproduction data in Ref. (3), but the un¬ 
certainties there were larger (above 40%). Our results thus confirm previous 
assertions that the LD corrections are at the few percent level only [3,4] and 
further show that these contributions are well under control. The amplitude 
for b —» sq including SD and LD contributions can be expressed as [3] 

A TOT (b -»• 57 ) = fV tb V t * s ^ m b Cf{m b ) - a 2 (m b )^— ^ 

• sa ltv RbF IJV . (11) 

where m s , M s have been neglected compared to m b . Using a- 2 (;m b ) ~ 0.24 ± 
0.04[3], Cj S (mb) = —0.30 ± 0.03 [16] and m b = 4.8 ± 0.2GeV, we find that 
the LD contribution increases the magnitude of the amplitude by (4 ± 2)%. 
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6. Application to B s - B\ 7 

Another process where the s —> dy quark transition will dominate is 

—> B d 7 . There are no pole contributions and we explicitly estimated the 
LD contribution from light meson loops to be smaller but comparable with 
the SD s —> dr/ contributions. This is an unusual B s meson decay in the 
sense that it represents the decay of the light quark in a Qq system. Also, it 
has a clear signature: two photons with energies of about 50MeV and 46MeV 
(the second one coming from the decay B* 0 \ —> Bd 7 ), followed by a usual B d 
decay. 

We roughly estimate the B s —> B d 7 decay rate from our s —> dy ampli¬ 
tude ( 6 ) by assuming that the spatial wavefnnctions of the s quark in the 
B s meson and the d quark in the B* d meson are similar, and noting that 
the photon energy (= 50MeV) is small compared to the average momentum 
(O(700MeV)) of the light quark in the bound state. A “free quark” approx¬ 
imation should then give a reasonable estimate of the transition amplitude. 
In terms of the effective s —» dq Hamiltonian ( 6 ) we obtain for the decay 
rate: 


r(B, B> dl ) 


a 


167T 4 


G 2 F \q\*{ m s F 2 + 


8n 2 va 2 CvMDMs 
Ml - M 2 


m d F 2 


87i 2 va 2 CyMDMd\ 
Ml - M 2 J 


( 12 ) 


where q is the photon momentum in the B s rest frame. Comparing to eq. (7) 
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and using the upper bound (1) we obtain 


T{B S -> B* 7) < 1.4 x 10 - 20 GeV . (13) 

Then, the present central value for the B s lifetime Tb s — 1-34 x 10 _12 s [14] 
gives a bound on the branching ratio, BR(B S —>• B d 7) < 3 x 1CT 8 . Although 
this is a very rare decay mode, its unique signature and the large number of 
B s mesons expected at B meson factories and at LHC-B 0(2 x 10~ n ) [23] 
make it interesting. 

7. Conclusions 

Using a VMD approximation, we found that the LD contribution to the 
s —> d'y transition may be significantly larger than the SD one, and could 
even lead to a saturation of the present experimental upper limit on the 
decay rate for fC —> (eq. (1)). This result throws new light on this 
decay mode. A further tightening of this upper limit or a measurement of 
the fl - —> rate would provide us with very useful information about 

the relative importance of the LD and SD contributions to s —> d'y. The 
present upper bound already implies a non-trivial cancellation at a level of 
30% or better in the LD contribution. The resulting approximate relation 
(eq. (3)) allowed us to estimate the relative importance of the LD contri¬ 
bution to the b —> 57 transition amplitude. Our estimate of (4 ± 2)% for 
this relative LD contribution agrees with earlier ones, which had larger un¬ 
certainties. Because the unusual process B s —>• B* d 7 is also dominated by an 
s —> d'y transition, its decay rate is related to that of O - —> S“ 7 . We find 
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that a present limit on the latter (eq. (1)) implies an upper bound for the 
branching ratio BR(B S —> B*^) < 3 x 10~ 8 , which is small but hopefully 
accessible in future experiments. 
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